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We analyze the low-lying spectrum and eigenmodes of lattice Dirac operators with a twisted mass 
term. The twist term expels the eigenvalues from a strip in the complex plane and all eigenmodes 
obtain a non-vanishing matrix element with 75. For a twisted Ginsparg- Wilson operator the spec- 
trum is located on two arcs in the complex plane. Modes due to non-trivial topological charge of 
the underlying gauge field have their eigenvalues at the edges of these arcs and obey a remnant 
index theorem. For configurations in the confined phase we find that the twist mainly affects the 
zero modes, while the bulk of the spectrum is essentially unchanged. 
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In the last few years lattice QCD formulations with a twisted mass term were studied in great 
detail |Jl]l-[|5|] and the latest results were reviewed at this conference [§]. The interest of the lattice 
community in twisted mass QCD is based on the fact that the twisted mass term protects the Dirac 
operator from small eigenvalues, thus avoiding exceptional configurations. This allows one to work 
at small quark masses and the results can be extrapolated to vanishing twist. 

However, the fact that no small eigenvalues appear implies that also zero modes of the Dirac 
operator are excluded. Zero modes of the Dirac operator in turn are, via the index theorem [|7|], 
related to infrared structures of the gauge field carrying topological charge. This connection already 
raises the question we address in this contribution: How does topological charge of the gauge field 
manifest itself in the eigensystem of the Dirac operator with a twisted mass term? 

We stress that this question is not only of academic interest, since topological objects are be- 
lieved to play an important role in the QCD vacuum. Thus it should be understood how topological 
objects interact with the fermionic degrees of freedom in the presence of a twisted mass term. 

For two mass-degenerate flavors of fermions, u and d, the twisted mass Dirac operator 9 has 
the form 

9 = (Do + m)t 2 + iMfi}T3 , t 3 = diag(l,-l), 

where I2 and T3 act in flavor space. For our analysis it is convenient to analyze the Dirac operators 
for u and d separately by setting D u = D(ju) and =D(— ju) with 

D(jU) = Do + m + inY 5 . (1) 

The Dirac operator Do is the 1 -flavor Dirac operator, discretized on the lattice without doublers. A 
common choice is Wilson's Dirac operator, but here we will also consider (for theoretical purposes) 
a Dirac operator that is chiral, i.e., obeys the Ginsparg -Wilson relation 

A property which we assume for the operator Do, and which is obeyed by Wilson's Dirac 
operator and also the overlap operator, is 75-hermiticity, 

7 5 Do75=Dj. (2) 

For the twisted operator D(/i) this property implies the generalized 75-hermiticity relation []|] 

7 5 DGu) 75 = D(- A t) t . (3) 

The generalized 75-hermiticity relation can be explored to derive results for the eigensystem of the 
Dirac operator. We denote an eigenvector of D(/j.) with eigenvalue A by yfy and find for its matrix 
element with 75 : 



KYx,JsWx) = (Vx,75D(n)Yx) = {VxM-V-YYsVx) = (D(-n)Yx,75Yx) = 
([D(n)-i2nY 5 ]Yx,Y5¥x) = ^*(Y}i,r5Yx)+i^(Yx,¥l) = ^*(¥x,Y5Yx)+i^- (4) 

We have used (|3|) and the fact that the eigenvectors are normalized to 1. From the last equation 
follows 

(VA,75VA) = ^. (5) 
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Figure 1: Left-hand side plot: Spectra of the Dirac operator in the complex plane for jj. = 0.0 and jl = 
0.02. The numbers next to the symbols for the eigenvalues are the matrix elements of the corresponding 
eigenvectors with 75. Right-hand-side plot: 75 matrix elements of the eigenvectors y/^ as a function of ImA. 
The full curves are the hyperbolas of Eq. (5). 



This formula implies that all eigenvectors have non-vanishing chirality which decreases monoton- 
ically with ImA. This is very different from the case without twist term, where only eigenmodes 
with real eigenvalue have a non- vanishing matrix element with 75 . These modes correspond to zero 
modes in the continuum and are of topological nature. 

Using the fact that 75 is bounded, i.e., \(YXiYs¥?i)\ < h we find (/I > 0) 

I ImA I > n . (6) 

Thus all eigenvalues have an imaginary part larger or equal the twist parameter jx (positive in our 
convention). The twist term expels the spectrum from a strip of width 2\i along the real axis and 
thus the smallest eigenvalues have a modulus with a size of at least \i. 

The basic properties of the spectrum and the eigenmodes are illustrated in Figure 1 which 
was generated using the twisted Wilson Dirac operator on a quenched gauge field configuration of 
topological charge Q = +1 on a lattice of size 16 4 at a lattice spacing of a = 0.094 fm as determined 
from the Sommer parameter JT0[]. 

In the left-hand side plot we compare the spectra in the complex plane for vanishing twist and 
at a twist of /I = 0.02 in lattice units. As the twist term is turned on, the eigenvalues are shifted away 
from the real axis. The topological eigenvalue (filled symbol) is shifted in the negative imaginary 
direction. It is obvious that the spectrum is no longer symmetric with respect to reflection at the 
real axis. We remark that the fluctuations of the eigenvalues in horizontal direction are not reduced 
by the twist term. The numbers next to the symbols representing the eigenvalues A are the values 
for the matrix elements (Y^,YsYx) of the corresponding eigenvectors Yx- 

In the right-hand side plot of Figure 1 we show these matrix elements as a function of ImA 
for three different values of fi. The symbols represent the numerical data, the full curves show the 
hyperbolas of Eq. (5). 
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As discussed, the focus of this paper is on the low lying spectrum of the twisted mass lattice 
Dirac operator and the role of the index theorem. For the index theorem it is important to use a 
setting where chiral symmetry is properly implemented on the lattice and to consider the role of 
the twisted mass term in this case. Thus we require now the Dirac operator Do (before the twist 
term is turned on) to obey the Ginsparg-Wilson relation, 

y 5 D + D 75 = aDoYsDo , (7) 

which governs chiral symmetry on the lattice. In this case D can have exact zero modes Yq which 
are chiral, i.e., JsYo = ^Vo- They are related to the topological charge Q of the gauge field via 
the index theorem, 

Q = n_ - n+ , (8) 

where n + and ra_ are the numbers of right- and left-handed zero modes JTT|]. 

Since the zero modes Yq °f A) are chiral, they are also eigenmodes of D(/i) with eigenvalue 
±ifl , with a plus sign for the right-handed mode and a minus sign for a left-handed zero mode. 
Thus in the presence of the twisted mass term we can identify a remnant of the index theorem 
which still has the form of Eq. (S), but and n + are now interpreted as the number of eigenvalues 
+i{JL and — ijX. 

From ([[]) one finds Dq = D(/j.) — //175 and inserting this expression into ^ one obtains after a 
few lines of algebra (use (||) and D(—\i) = D(/i) — i2jj,y 5 ), 

D(jii) + D(/i) t = aD(v)D{^ ~ a\l 2 . (9) 

Sandwiching this equation between eigenvectors Yk gives rise to the equation 

X +X* = aXX* -a/x 2 , (10) 
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Q = -1 , real Polyakov loop 



Q = +1 , complex Polyakov loop 
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Figure 3: Spectrum of the Dirac operator in the deconfined, chirally symmetric, phase. We compare the 
untwisted case to the situation at jj. = 0.02. The left-hand side plot is for a configuration with Q = — 1 and 
real Polyakov loop, while the right-hand side plot is for Q = 1 and complex Polyakov loop. 



for the corresponding eigenvalue A. Setting X = x + iy one finds that this is the equation for a 
circle in the complex plane. This "stretched Ginsparg-Wilson circle" has center 1 /a and radius 
y/ 1 1 a 1 + /I 2 . Since \Imk \ > /I we find that the spectrum of D(/i) is confined to the two arcs shown 
in the left-hand side plot of Fig. 2. The topological eigenmodes with eigenvalues ±ZjU fall one the 
edges of the two arcs (on the other edges their doubler partners are found). The matrix elements of 
the eigenvectors with 75 still obey Eq. (|5]). The right-hand side plot of Fig. 2 illustrates the results 
for the 75 -matrix element. 

Our discussion of a Ginsparg-Wilson operator with a twisted mass term shows that the zero 
modes of a chiral Dirac operator without twist, Do, are protected by chiral symmetry and remain 
unchanged as the twist term is turned on (only the eigenvalue changes). Instead, when using the 
Wilson Dirac operator for D , this protection is only approximate. In [^] characteristic properties 
of the topological zero mode (edge mode), such as its localization and 75 matrix element, were 
studied for different values of }X. While for the twisted Wilson operator these properties changed 
considerably, for the approximately chiral CI operator [12] they remained essentially unchanged. 

Another question that was addressed in []|] is the role of the twisted mass term in the de- 
confined (chirally symmetric) phase. In the chirally symmetric phase the spectrum of the Dirac 
operator has a gap near the origin (up to isolated zero modes) which implies, via the Banks-Casher 
relation [13], a vanishing chiral condensate. Turning on the twisted mass term mainly affects the 
zero modes. As seen in Fig. 3, they are shifted up or down, depending on their chirality. The size 
of the spectral gap remains essentially unchanged. 
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